Genuine Multipartite Entanglement in Quantum Phase Transitions 
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We demonstrate that the Global Entanglement (GE) measure defined by Meyer and Wallach, J. 
Math. Phys. 43, 4273 (2002), is maximal at the critical point for the Ising chain in a transverse 
magnetic field. Our analysis is based on the equivalence of GE to the averaged linear entropy, 
allowing the understanding of multipartite entanglement (ME) features through a generalization 
of GE for bipartite blocks of qubits. Moreover, in contrast to GE, the proposed ME measure 
can distinguish three paradigmatic entangled states: GHZm, Wn, and EPR® N / 2 . As such the 
generalized measure can detect genuine ME and is maximal at the critical point. 
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Entanglement is a correlation of exclusively quantum 
nature present (in principle) in any set of post-interacting 
quantum systems [l( . As such multipartite entanglement 
(ME) is expected to play a key role on quantum phase 
transition (QPT) phenomena in the same way that (sta- 
tistical ^classical correlation does on classical phase tran- 
sitions [si Id]- In ordinary phase transitions, at the critical 
point, a non-zero order parameter characterizes a long 
range correlation (given by the correlation length diver- 
gence). In the same way, in QPTs it is expected that 
ME be maximal at the critical point, in the sense that 
all the system parties would be entangled to each other 
0. However, this conjecture could not be proved in gen- 
eral neither by measures of pairwise entanglement nor by 
the proposed ME measures. Even after a considerable 
effort, a deep understanding of multipartite entangled 
states (MES) is lacked. It is still a great challenge thus 
to capture the essential features of genuine ME, from a 
conceptual point of view, as well as from a quantitative 
approach, defining a measure that among other proper- 
ties be able to distinguish MES 

Indeed, concerning the legitimate quantum correla- 
tions in QPTs it would be certainly important to know 
exactly what kind of entanglement should we expect to be 
maximal at the critical point. The great majority of ef- 
forts trying to answer this question made use of two kinds 
of bipartite entanglement measures, both calculated for 
spin-1 /2 lattice models such as the Ising model in a trans- 
verse magnetic field Q|. The first one, namely the pair- 
wise entanglement (concurrence) between two spins in 
the chain, was studied by Refs. [^, 0|. The second 
one, the entropy of entanglement between one part of 
the chain (a block of L s_Dins)and the rest of the chain, 
was investigate by Refs. (a,tD,lq|- Some candidates of ME 
measures were also evaluated in systems exhibiting QPTs 
ISEHEIl Nevertheless, none of the entanglement mea- 
sures employed in the above references are maximal at 
the critical point but the single site entropy for the Ising 
model [a| in the thermodynamical limit and the Locafiz- 
able Entanglement 0] for an Ising chain with 14 spins. 
Furthermore, in Refs. 0, 0] the authors have indepen- 



dently shown that bipartite entanglement vanishes when 
the distance between the two spins is greater than two 
lattice sites 0|. This is not expected since long range 
quantum correlations should be present at the critical 
point. It was then suggested that bipartite entanglement 
at the critical point would be decreased in order to in- 
crease ME due to entanglement sharing jjj. In other 
words, ME only appears at the expense of pairwise en- 
tanglement and at the critical point we should expect a 
genuine MES. 

In this paper we demonstrate that the Global Entan- 
glement (GE) introduced in Ref. [13] indeed captures the 
essential point to be maximal at the critical point for the 
Ising model in a transverse magnetic field in the ther- 
modynamical limit. We also prove that there exists an 
interesting relation among GE, von Neumann entropy, 
linear entropy (LE), and 2-tangle [lH. flfil Hflj. showing 
that they are all equivalent to detect QPTs. Furthermore, 
this relation helps us to understand the results obtained 
in Ref. 0, as outlined in the previous paragraph, and 
suggests that they are not particular to the Ising model 
but common to all MES with translational invariance. In 
addition to this, we generalize GE and propose a new ME 
measure, which is also maximal at the critical point for 
the Ising model, can detect genuine MES, and contrary 
to GE, furnishes different values for the entanglement of 
the GHZ N , W n , and EPR® N / 2 states, thus being able 
to distinguish among MES. 

For a N qubit system (spin- 1/2 chain) it was noticed 
that GE is simply related to the TV single qubit purities 
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S L ( Pj ) = (S L ), (1) 



where GE is here on identified as Eq\ pj = Trj{p} is 
the j-th qubit reduced density matrix obtained by tracing 
out the other j qubits, and Sl(Pj) = g4i [l — Tr (pf)] is 
the standard definition of LE. This relation shows that 
Eq^ is just the mean of LE. It was also noticed in Refs. 
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J=l 

where rarest = G 2 is the 2-tangle [3, EH Oil- the 
square of the concurrence G Both LE and the 2- 

tangle can thus be used to quantify the entanglement 
between any block bipartition of a system of N-qubits. 
(They quantify the entanglement between one qubit j 
and the rest N — 1 qubits of the chain 0.) The proof 
of l|2j is based on the Schmidt decomposition |2lJ , which 
also allows us to use for pure systems the reduced von 
Neumann entropy, S v (pjQ)) = -Tr^j) [p Q) togApjQ))] , 
as a good bipartite entanglement measure [22j. Here 
d = min{dim7ij, dimHj} and dimTlj(j' ) is the Hilbert 
space dimension of subsystem j(j). Recalling that SV is 
bounded from below by Sl and employing Eqs. Q and 
(j2J) we obtain the following important relation 

4 X) = M = ( S l) < ( S v), (3) 

which states that GE is nothing but the mean LE of sin- 
gle qubits with the rest of the chain. Furthermore, GE 
is also equal to the mean 2-tangle and a lower bound for 
the mean von Neumann entropy. An immediate conse- 
quence of this result shows up when we deal with linear 
chains with translational invariance. This implies that 
(Sl) = Sl(Pj) and that (SV) — Sv(pj)- Hence, Eq. <EJ 
becomes E G — Sl{Pj) < Sv(pj)- Since Sl(pj) and 
Sv(pj) have the same concavity and both entropies at- 
tain their maximal value for a maximally mixed state this 
last relation shows that E G ^ is as efficient as the linear 
and the von Neumann entropies to detect QPTs. In Ref. 

the authors used Sv and in Ref. E G ^ was em- 
ployed to detect QPTs in the Ising model. Needless to 
say, both works arrived at the same results for a given 
range of parameters via, notwithstanding, different en- 
tanglement measures which by that time were thought 
to be unrelated. 

Despite its success to detect the Greenberger-Horne- 
Zeilinger (GHZ) state [HJ E^, E { G ] sometimes fails 
for distinguishing different multipartite states. This is 
best understood if we study Eq for three paradig- 
matic multipartite states. The first is \GHZn) = 
(l/V2)(|0)® JV + |l)® Ar ), where |0)®* and rep- 
resent N tensor products of |0) and |1) respectively. 
The second is a tensor product of N/2 Bell states [18| . 
\EPR N ) = |$+)®tt, where |*+) = (1/V2)(|00) + |11)). 
This state is obviously not a MES. Only the pairs of 
qubits (2j — 1,2 j), where j — 1,2, ...,N, are entangled. 
Nevertheless, for both states E G ^ = 1. The last one is 
the W state 0: \W N ) = (1/y/N) J2f=i |00- • • lj • • • 00). 
The state |00 • • • lj ■ ■ ■ 00) represents N qubits in which 
the j-th is |1) and the others are |0). As shown in Ref. 
0, E^(W N )=A(N -1)/N\ 



We now present a generalization of GE. The main fea- 
tures of this new approach are three-fold. First, it be- 
comes clear that we have different classes of ME mea- 
sures, where Eq is the first one. Second, the first non 

(2) 

trivial class, E G , furnishes different values for the three 
states considered above Third, it gives new insights in 
the study of QPT and ME. 

In order to define E G we need the following function 

A ( 1 N ~ 1 \ 

gm = 3 (i-^E^L+Oh (4) 

where Pjj+i is the density matrix of qubits j and j + I, 
obtained by tracing out the other N— 2 qubits. The index 
< I < N is the distance in the chain of two qubits and 
4/3 is a normalization constant assuring G(2,l) < 1. Of 
interest here are two quantities that can be considered 
ME measures in the same sense that E G ^ is: 

/ 1 N- l \ 

G(2,l) - j l-]^E*i) - ( 5 ) 

and jv_! 

4 2) = ^EG(2,0. (6) 
(=1 

We can interpret G(2, 1) as the mean LE of all two qubit 
nearest neighbors with the rest of the chain. Similar in- 

(2) 

terpretations are valid for the others G(2,l). E G is the 
mean of all G(2, 1) and it gives the mean LE of all two 
qubits, independent of their distance, with the rest of the 
chain To define E G ^ we need the function G(3, 1%, 1%) 
with one more parameter, since now we can have different 
distances among the three qubits of the reduced state. A 
complete analysis of this new ME measure and its use- 
fulness to detect MES is discussed elsewhere |2^ |. 

Table []] shows the quantities given by Eqs. JU and 
© for GHZ N ,EPR N , and W N . We note that due to 
translational symmetry, G(2, 1) and E G ' are identical for 
GHZn and Wn. It is worthy of mention that depending 
on the value of N, the states are differently classified by 
G(2,l). A similar behavior is observed for E G ^ [24j]. In 
this case, however, EPRn is the most entangled state for 
long chains. The reason for that lies on the definition of 
E { G \ For EPR N , G(2,l) = 1 for any I > 2. Thus, since 
Eq is the average of all G(2, 1), for long chains G(2, 1) 

(2) 

does not contribute significantly and E G -» 1. 

It is worth noticing that even at the thermodynamical 

( 2) 

limit, N — > oo, E G and G(2, 1) still distinguish the three 
states. However, the ordering of the states is different. 
As already explained, this is due to the contribution of 
G(2,0, I > 2, in the calculation of E { G ] (EPR N ). 

Now we specify to the one-dimensional Ising model in a 
transverse magnetic field, which is given by the following 
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Table I: Comparison among the three paradigmatic states. 
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for A > 1. Here 



l u l I — \"1 

for A < 1 or 

tt JO Ufb l + A 2 +2Acos(fc) ' 

and I > 1 is the lattice site distance between two qubits. 
By tracing out one of the qubits we obtain the single 
qubit density matrix, which allows us to obtain E^ as 
a function of A. This is shown in Fig. As a matter 
of fact E^q is maximal (with singular derivative) at the 
critical point A = 1. For comparison, in Fig. 1 we plot 
Sv(pj), which was already shown also maximal at the 
critical point for the broken-symmetry state We em- 



where i represents the i-th qubit, A is a free parameter re- 
lated to the inverse strength of the magnetic field, and we 
work in the thermodynamical limit. We assume periodic 
boundary conditions: <7/v+i = a\. As we have shown, for 
a system with translational symmetry GE is nothing but 
LE of one spin with the rest of the chain. We only need, 
then, LE to obtain GE. For that end we must calculate 
the single qubit (or single site) reduced density matrix, 
which is obtained from the two qubits (two sites) reduced 
density matrix. It is a 4 x 4 matrix and can be written 
as 

Pij = r ^ij\p] = \ J2 Pa l 3<T ? ® °f ' 

where p is the broken-symmetry ground state in the ther 
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modynamical limit and p a p 



^[afa p jPij ] 



(aro 



0\ 



Tijj is the partial trace over all degrees of freedom except 
the spins at sites i and j, af is the Pauli matrix acting on 
the site i, a, (3 — 0, x, y, z where a is the identity matrix, 
and Pap is real. Therefore, all we need are the ground 
state two-point correlation functions (CFs). By symme- 
try arguments concerning the ground state Q the only 
non-zero CFs are p 00 , p xx , P yy , Pzz, Pox = PxO, Poz = Pzo, 
and pxz — Pzx- Due to normalization poo = 1 and a 
direct calculation gives p xz = p zx = for A < 1. On 
the other hand, the Schwartz inequality necessarily gives 
< \p X z\ < \( (J i)(< J i)\, allowing thus that the lower and 
upper bounds for entanglement be calculated for A > 1. 
We plot the upper bound for entanglement by taking 
Pxz = 0. By continuity the true value for entanglement 
must show a similar behavior. 

Those CFs have been already calculated 0| and we just 
highlight the main results. The two-point CFs and the 
mean values of a x and a z are 
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Figure 1: (Color online) Von Neumann entropy (dashed) and 
GE/LE (solid) as a function of A. 

phasize that these measures quantify entanglement in the 
global system by measuring how mixed the subsystems 
are. The physical meaning behind studying "mixedness" 
lies on the fact that the more entangled two subsystems 
are the more mixed their reduced density matrix should 
be 0, 0| . However, in a many-body system there are 
many ways in which one could divide the global system 
into subsystems. The first non-trivial generalization is 
to study LE of two sites with the rest of the chain. Us- 
ing we can calculate G(2, 1) for the Ising model (Fig. 

El- It has a similar behavior to Eq\ being also maximal 
(with singular derivative) at the critical point. This fea- 
ture demonstrates that both a pair of nearest neighbors 
sites and the sites themselves are maximally entangled to 
the rest of the chain at the critical point. But this is not 
particular to nearest neighbors as shown in Fig. UH where 
G(2,l), G(2,15), and E^ ] = ^ Ei=i G(2, i) is plotted. 
G(2, 15) is also maximal at the critical point, indicat- 
ing that in a QPT entanglement sharing at the critical 
point is favored by an increase of all kind of ME. More- 
over, Fig. El shows that G(2, 15) is only slightly different 
from Eq^ = A Yll=i G(2, i) . This is due to the rapid 
convergence of G(2, 1) as I is increased. At the critical 
point lim/^oo G(2, 1) is 0.675, and thus higher than the 
value for GHZn, EPRn, and Wn, obtained in the ther- 
modynamical limit, indicating thus a genuine MES. We 



also remark that besides Eq \ G(2, 1), and Eq' being all 
maximal at the critical point, E^ < E^ for every value 
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of A. However an interesting change of ordering for E G 
and G(2, 1) occurs around the critical point. For A < 1, 
E { q ] > G(2,l), but for A > 1, < G(2,l). Thence 
a kind of ME is favored in detriment of the other, de- 
pending on the system phase. Also, the fact that at the 
critical point both E^ and Eq are maximal indicates 
entanglement sharing, such that all the sites of the chain 
are strongly (quantum) correlated. Of course this state- 
ment is only completely true provided that Efc' is also 
shown to be maximal for any 2 < m < N — 1 (all possi- 
ble partitions). Furthermore, the fact that G(2, 1) always 
increase as I — > oo at the critical point suggests a kind 
of diverging entanglement length. However its precise 
definition demands a careful calculation of the scaling of 
entanglement such as in Refs. 00- These points are 
left for further investigation [25]. 
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Figure 2: (Color online) E*q (solid) and G(2, 1) (dashed) as 
a function of A. Both quantities are maximal at the critical 
point A = 1. 
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Figure 3: (Color online) G(2, 1) (dashed/black), G(2, 15) 
(solid/red), and Eq (dotted-dashed/blue) as a function of 
A. We see that E G is slightly different from G{2, 15), show- 
ing that G(2, 1) saturates as I — » oo. 

In conclusion we have demonstrated that for an infi- 
nite Ising chain both Eq and its generalization, Eq\ 

(2) 

are maximal at the critical point. Furthermore, E G ' as 
defined here is able to detect genuine ME. We remark 
that the behavior of the ME measures here presented for 
an infinite chain is in agreement with the Localizable En- 



tanglement calculated for a finite (N=14) Ising chain for 
the broken-symmetry state ^lj • Yet our results were ob- 
tained in a relatively simpler fashion and could be used 
to infer genuine ME for systems where the Localizable 
Entanglement has failed to detect QPT_|2jj] . Finally, our 
results reinforced Osborne and Nielsen [5J conjecture that 
at the critical point ME should be high, due to entangle- 
ment sharing, in detriment of bipartite entanglement. 
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